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Efimov physics in two nuclear-spin sublevels of bosonic lithium is studied and it is shown 
that the positions and widths of recombination minima and Efimov resonances are identical 
for both states within the experimental errors which indicates that the short-range physics 
is nuclear-spin independent. We also find that the Efimov features are universally related 
across Feshbach resonances. These results crucially depend on careful mapping between 
the scattering length and the applied magnetic field which we achieve by characterization 
of the two broad Feshbach resonances in the different states by means of rf-spectroscopy of 
weakly bound molecules. By fitting the binding energies numerically with a coupled channels 
calculation we precisely determine the absolute positions of the Feshbach resonances and the 
values of the singlet and triplet scattering lengths. 



I. INTRODUCTION 

The recent years' remarkable progress in the study of Efimov quantum states in ultracold 
atoms has renewed a great deal of interest in this "old- new" quantum few-body problem. Since 
the first prediction made by V. Efimov in the early 70's [1], many systems were considered for 
an experimental study of these quantum states, however all of them were attempted in vain. The 
first experimental evidence of Efimov physics was reported in 2006 in a system of ultracold ^'^^Cs 
atoms which was later on enhanced and verified in an additional study of the same system . 
Since then, signatures of Efimov physics have been observed in other ultracold atomic species which 
turn out to be the only platform up to now suitable to study Efimov physics. 

It is well known that at very low collision energies the only partial wave contributing to the 
scattering process is the s-wave 

and thus the two-body interaction is completely determined by the s-wave scattering length a. 
When a exceeds the characteristic two-body potential range tq, weakly bound three-body Efimov 
states emerge and their number scales logarithmically with a, N}j = (sq/vt) ln(|a|/ro) where sq = 
1.00624 m. When I a I — )• iboo, the number of bound states goes to infinity and their energies are 

n R 

related in powers of a universal scaling factor exp(— 27r/so) ~ 1/(22.7)^ [ll, Uj- A first indirect 



evidence of two consecutive Efimov states was demonstrated in ultracold ^^K [al and a universal 

n 

scaling across a region of \a\ — t- ±00 was verified in "^Li [6]. 

Studies of Efimov physics have been rapidly extended beyond the spin polarized bosonic samples. 
A notable example of a three-fermion system, all in different spin states, has been shown by a 
number of experimental groups 7h9(] one of which developed a new and promising experimental 
app^ch to p„be d.ectly the Efi.ov ,™„tu,. st.t. Q. Another example of a d.ffetent system 
is heteronuclear universal trimers observed in Ref. 

Recent developments extended the universal few-body physics to the domain of four-body states 



which was theoretically predicted 
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and experimentally verified jisl ] . Many three- and four-body 



features over a large dynamical range have been demonstrated in ^Li in Ref. [l^- Interestingly, 
the universal few-body physics can .ee.te„«toN..oC.c,tt.te.a.hott.ht;:.e.pe— 
confirmation remains obscure. 

In this work we summarize our study of Efimov physics in two different energy sublevels and 
across two different Feshbach resonances of the same atomic system P, [l^- As both Feshbach 
resonances occur at high magnetic fields where nuclear and electronic spins are decoupled, the two 
energy sublevels are associated with two nuclear-spin states. Our main finding is that the Efimov 
features are identical within the experimental errors in these two states. The absolute location and 
lifetime of an Efimov state is defined by the unknown short-range part of the three-body potential. 
Most generally, the short range potential is given in terms of two-body potential permutations 
of the two-body subsystems and a true three-body potential which is of importance only when 
three particles are very close together. It is very difficult to solve the short-range physics exactly, 
and therefore this region is usually treated in terms of a three-body parameter J, 17]. Thus, our 
results should be interpreted as a proof that the three-body parameter (and thus the short range 
physics) is identical for the two states. We provide new and accurate characterizations of Feshbach 
resonances on both states and reevaluate our previously published results in accordance with this 
study. Though small changes in the positions of the Efimov features can be identified, the main 
conclusion is not affected. 

This paper is organized as follows. Study of Efimov physics across two different Feshbach res- 
onances in two nuclear-spin states is reported in Section |TI1 We start by comparing the relevant 
properties of the two energy levels and their resonances. Then, after a description of the experimen- 
tal procedure we report on newly fitted positions and widths of the Efimov features based on our 
most recent and precise characterization of the Feshbach resonances on both nuclear-spin states. 
The latter is based on fitting of the binding energies of weakly bound molecular states , obtained by 
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RF spectroscopy, to a coupled channels calculation locally for each resonance and globally for both 
resonances. This study is crucial for an accurate mapping of the magnetic field to the scattering 
length and it is considered in Section HIH Finally, Section [TVl concludes the paper. 



II. EFIMOV PHYSICS ACROSS FESHBACH RESONANCES 
A. Comparison between two nuclear-spin states 




Figure 1: (a) The broad Feshbach resonances of the \mp = 1) (dashed red line) and the \mp = 0) (solid 
blue line) states, centered at 738.2 G and 893.7 G, respectively, (b) The two Feshbach resonances of the 
Imp = 0) state. 

A comprehensive discussion on the experimental characterization of Feshbach resonances will 
be the subject of Section IIIII Here we intend to draw a general comparison between the broad 
resonances of the Imp = 1) and the Imp = 0) states of ''Li in the context relevant to the study of 
Efimov physics. In Fig. [DJa) we show the scattering length in units of Bohr radius gq as a function 
of magnetic field in the vicinity of both Feshbach resonances with their centers aligned. It can be 
easily recognized that the resonances are comparable in widths while the one on the Imp = 0) 
state is slightly wider. Although there are different hyperfine states involved, these two resonances 
have their origin in the same molecular bound state. Another difference is that in contrast to the 
Imp = 1) state where only one resonance exists, on the \mp = 0) state there is a second resonance 
overlapping with the first one. It is much narrower and is positioned by atom loss measurement at 
~ 845 G [a] as shown in Fig. [U^b). We note that the measurements reported here are obtained in 
close vicinity to the wide resonance's position, away from the narrow one by many times its width. 
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thus it is not expected to influence the resuhs in the region of interest. 

Signature of Efimov physics is studied here by measuring three-body recombination loss of atoms 
as it has been investigated in all but one recent experiments. In that respect, there is an inherent 
difference between the two states: in the absolute ground state, two-body relaxation mechanism is 
fundamentally forbidden while in the {mp = 0) state it is allowed. In Fig. [2] we show the dipolar 
relaxation rate coefficients as a function of magnetic field which were calculated via a coupled- 

n 

channels calculation by using recent interaction potentials [18]. Except for two peaks which signify 
Feshbach resonances, the loss rate coefficients are extremely small, ~ 3 orders of magnitude smaller 
than the corresponding measured rate coefficients, if the experimental losses were treated as purely 
two-body related. For example, at 880 G the 2-body loss coefficient is ~ 5 * 10^^^ cm^/s which 
yields, given that the atom density is ~ 10^^ cm~^, a life time of ~ 20000 s. As a comparison, the 
life time in our dipole trap due to vacuum is less than 100 s. As a result, we exclude two-body losses 
from the analysis described below and determine that the loss processes in the region of interest 
are related to three-body recombination. Note that while we find this mechanism to be negligible 
for ^Li atoms, it can be important for heavier alkali atoms. For instance, ^^^Cs experiences large 
dipolar losses caused by the second-order spin-orbit interaction [l^. 
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Figure 2: A coupled-channels calculation of dipolar relaxation rate coefficients as a function of magnetic field 
for the Imp = 0) state. 
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B. Experimental measurement of three-body recombination loss 

Three-body recombination loss rate near a Feshbach resonance has a general dependence and 
in the zero-temperature limit it diverges at the resonance [20|] . Above this general scaling, universal 



theory predicts log-periodic oscillations of the three-body recombination loss rate coefficient {K3) 
due to the presence of Efimov trimer states. For positive scattering lengths the oscillations are 
caused by destructive interference conditions between two possible decay pathways at certain values 



of a 
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2l| . For negative scattering lengths the loss rate coefficient exhibits a resonance enhancement 



each time an Efimov trimer state intersects with the continuum threshold. Hence, a study of 
across a Feshbach resonance enables the search for an indirect evidence of Efimov trimer states. 



25000 



I 20000 
o 

CI) 

g 

15000 



E 
0) 



10000 




0.1 0.2 0.3 0.4 0.5 
♦ ♦ Time (msec) 




0.05 0.10 
Time (sec) 



Figure 3: A typical atom-number decay measurement on the \mp = 1) state. Inset: time-of-flight measure- 
ment from which the temperature is deduced. 



In the experiment, we perform evaporative cooling in an optical dipole trap near a Feshbach 
resonance where a gas of ^Li atoms is spontaneously spin purified to the Imp = 0) state and is 
cooled down to the threshold of degeneracy [S]. For measurements in the absolute ground state 
{\mp = 1)), a rapid adiabatic passage by means of a radio- frequency (rf) sweep at a low magnetic 
field is used to transfer the atoms from the \mp = 0) state jisl . We conduct measurements of 
atom-number decay and initial temperature as a function of magnetic field in the vicinity of wide 
Feshbach resonances on both states. An example of such set of measurements is shown in Fig. [3l We 
use the temperature measurement, together with a precise characterization of the trap frequencies, 
to estimate the initial atom density. Then, the value is extracted by fitting the atom-number 
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decay measurement with the atom-loss rate equation solution: 

N = -K3{n'^)N -TN, (1) 

where and T are the three- and single-body loss rate coefficients, respectively. T is determined 
independently by measuring a very long decay tail of a low density sample. Note that this simplified 
model does not include several important effects, one of which is the saturation of to a maximal 
value Kmax due to finite temperature (unitarity limit) and it can be represented as |23l |: 



where m is the atomic mass, /c^ is the Boltzmann's constant and c is a numerical constant distin- 
guishing between the two threshold regimes of the collision energies, c = 1 for a > and c ~ 0.1 
for a < 0. To address this limitation, measurements for which Kj, > O.lKmax are omitted from the 
analysis of Efimov features (see Section III Cp . 

Other effects which are not included in the model are 'anti-evaporation' and recombination 
heating 



24( 1 . The first is an effective heating caused by a preferential loss of atoms from the densest 
(and thus coldest) part of the cloud. We treat the evolution of our data to no more than 30% 
decrease in atom number, as can be seen in Fig. [21 for which 'anti-evaporation' is estimated to 
induce a systematic error of ~ 23% towards higher values of K^. The error is evaluated based on 
an analytical solution to coupled atom-loss rate and temperature evolution equations and it is 

estimated not to limit the accuracy of the reported results. Recombination heating can be neglected 
for K-^ <^ Kmax because in this case the trap depth (which scales with temperature) is smaller 
than the energy released in the recombination process leading to an immediate loss of the colliding 
partners. Moreover, this approach, together with the small decrease in atom number, allows time 
evolution of the unitarity limit (Kmaxit)) to be neglected within a single atom- number decay 
measurement. Note, that 'anti-evaporation' and recombination heating pose strong limitations on 
the measurement's dynamical range in the large scattering lengths limit where more careful model 
than Eq.Q has to be considered. This is demonstrated for \mp = 1) state in Fig. [4] where a 
dramatic increase in the initial temperature of the atom-number decay measurement near the 
resonance corresponds to a reduction in the value of Kmax according to Eq.Q. 

C. Efimov features and universality 

Experimental results of the three-body loss measurements are summarized in Fig. [5] where 
is plotted as a function of the scattering length a for the \mF = 1) state (red solid circles) and the 
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Figure 4: Initial temperature of the atom-number decay measurement around the Feshbach resonance on the 
\mF = 1) state from which the atom density is extracted for the calculation of K^. Each point is deduced 
from a time-of- flight measurement (see inset in Fig. [3l). The error bars represent the fitting errors. 



\mp = 0) state (blue open diamonds). Magnetic field values are converted into scattering lengths by 
using carefully characterized Feshbach resonances, discussed in Section IIIII A qualitative analysis 
indicates a striking similarity between the two sets of measurements for both positive and negative 
scattering lengths. We can further verify this similarity quantitatively by fitting the measured 
data to a prediction of universal theory. For that purpose we represent the loss rate coefficient in 
a convenient form [4]: 



= ?>C±{a)ha^/m 



(3) 



where it hints at the positive (+) or negative (-) region of the scattering length. An effective field 
theory provides analytic expressions for C±{a) that we use in the form represented in 



C+{a) = 67.1e-2'^+(cos2[soln(a/a+)] +sinh2r?+) + 16.8(1 - e 



(4) 



and 



C_(a) = 4590sinh(2r?_)/(sin2[soln(|a|/a_)] +sinh2ry_), 



(5) 



where the free parameters of the fits are a± and r]± which are connected to the real and the 



imaginary parts of the three-body parameter, respectively 



by solid lines in corresponding colors in Fig. [5] 



251 ] . The fitting results are represented 



26l | and the fitting parameters are summarized in 
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Figure 5: Experimentally measured three-body loss coeiBcient K3 as a function of scattering length (in units 
of Bohr radius ao) for the Imp = 1) state (red solid circles) and the Imp = 0) state (blue open diamonds). 
The solid lines (red and blue, respectively) represent fits to the analytical expressions of universal theory 
[2^. The error bars consist of two contributions: the uncertainty in temperature measurement which affects 
the estimated atom density and the fitting error of the atom- number decay measurement. 

Table HI Comparing corresponding parameters in different states acknowledges the notion that the 
Efimov features are identical within the experimental errors which signifies that the short-range 
physics is nuclear-spin independent. An interesting conclusion can be drawn from this observation. 
Recall that the short range potential is given in terms of two-body potential permutations of the 
two-body subsystems and a true three-body potential. As the two Feshbach resonances occur at 
high magnetic fields where the nuclear and electron spins are effectively decoupled, the two-body 
potentials are similar for both states. Therefore, if the short-range physics is spin- independent, the 
true three-body forces are either also nuclear-spin independent or play a relatively minor role. 



state 




77+ 77- a+/ao 


a-/ao 


a+/\a-\ 


Imp = 


0) 


0.213(79) 0.180(48) 238(25) 


-280(12) 


0.85(11) 


Imp = 


1) 


0.170(41) 0.253(62) 265(16) 


-274(12) 


0.97(8) 



Table I: Fitting parameters to universal theory obtained from the measured K3 values of the \mp ~ 1) and 
the Imp — 0) states. 

Moreover, two predictions of universal theory are also verified here: first, the decay parameters 
r/_|_ and r/_, which describe the lifetime of the Efimov state, are assumed to be equal and indeed 
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they are (within the experimental errors). This suggests that the imaginary part of the three-body 
parameter is identical. Second, the theoretical assumption that the real part of the three-body 
parameter across a Feshbach resonance is the same for negative and positive scattering length 
regions requires a+ and a_ to obey a universal ratio a+/|a_| = 0.96(3) [4]. The fits yield values of 
0.85(11) and 0.97(8) for {mp = 0) and Imp = 1), respectively, which overlap with each other and 
with the predicted value within the experimental and theoretical error bars. We thus confirm that 
the three-body parameter is preserved across two different Feshbach resonances and between two 
different nuclear-spin states. 

For positive scatteringlengths, the Efimov trimer is expected to intersect with the atom-dimer 
threshold at a* l.la+ J]. Theory predicts that a* and a_ of the same trimer state are related 



as a_ ~ —22a* Therefore, if the observed Efimov resonance at a_ indicates the lowest state, 
the one expected at a* indicates the first excited state as the lowest one becomes nonuniversal. 

III. CHARACTERIZATION OF FESHBACH RESONANCES 
A. Experimental measurements of binding energy of Feshbach molecules 

For a correct investigation of the coefficients dependence on the two-body scattering length, 
it is crucial to have an accurate mapping between the scattering length and the applied magnetic 
field. For this purpose, the binding energies of the underlying bound states which give rise to the two 
broad Feshbach resonances are carefully measured using rf molecule association. This method uses 
a weak rf field to resonantly associate weakly bound Feshbach dimers which are then rapidly lost 
through collisional relaxation into deeply bound states j^]- In the experiment the rf modulation 
time is varied between 0.5 and 3 sec and the modulation amplitude ranges from 150 to 750 mG. The 
remaining atom number is measured by absorption imaging as a function of rf frequency at a given 
magnetic field and the rf-induced losses are then numerically fitted to a convolution of a Maxwell- 
Boltzmann and a Gaussian distributions. The former accounts for broadening of th e sp ectroscopic 
feature due to finite kinetic energy of atoms at a typical temperature of ~ l.BfiK [28\. The latter 
reflects broadening due to magnetic field instability and shot-to-shot atom number fiuctuations. 
From the fit we extract the molecular binding energy (Ei,) corresponding to zero temperature. An 
example of a molecule association induced loss feature is depicted in Fig. [6] where the characteristic 
asymmetry of the obtained profile is clearly seen 23] . Results of the binding energy rf spectroscopy 



for both states are shown in Fig. [71 We analyze these results by fitting them with a coupled-channels 
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Figure 6: rf association of molecules at B — 734.4 G on the |m_F — 1) state. The loss resonance is fitted 
numerically to a convolution of Maxwell-Boltzmann and a Gaussian distributions (solid line) . 
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Figure 7: rf spectroscopy of the molecular binding energy near the Feshbach resonance in the \mp = 0) (a) 
and the \mF = 1) (b) states. The solid line represents an independent fit to CC calculation while the dark 
dashed line is the combined fit for both states simultaneously. 



B. Analysis of two-body interactions 



Very close to resonance, the binding energy is related to the scattering length as 



Ebia) 



(6) 
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This expression shows that ultracold scattering physics is closely related to bound state physics just 
below the dissociation threshold. However, this simple relationship rapidly breaks down when going 
further away from resonance, and is therefore not sufficient to interpret three-body recombination 
in terms of the scattering length. 

In order to make a correct mapping between binding energy and scattering length, we make use 
of a coupled-channels two-body interaction model for lithium, which has been discussed earlier [l^, 
and we use the interaction potentials discussed there as a starting point for our analysis. For large 
interatomic separations r the singlet (S = 0) and a triplet {S = 1) potentials are given by 



Vs{r) 



_ Cs_ _ Clo _ ( _-.xS f-iex 7/2a-l -2ar 
^6 ^8 ^10 ^ ^ ^ ' 



(7) 



3^, 



33] , and the ionization 



with van der Waals coefficients Cq,Cs, Cio [311, exchange parameter 
energy [34]. For the short radial range r we use model singlet and triplet potentials which 



have also been used in Refs. 



35|, 



The short-range and long-range potentials are smoothly connected at r = ISoq. We improve 
the accuracies of the short range potentials considerably, by making use of the accumulated phase 
method 35|. Therefore a boundary condition is applied on the partial- wave radial wave functions 



at r = 7ao in the form of a WKB phase 



(^s,t{E, £) = <PIt{E, £) + A(/>5,T. (8) 

The first term on the right is calculated by radial integration of the model potential up to 
7ao and is sufficiently accurate to account for the energy and angular momentum dependence of 
the accumulated phase. The corrections A(ps,T to the accumulated singlet and triplet phases are 
independent of energy and angular momentum. 

The most crucial parameters in the coupled channels model are A(f)s, and Cg, which we 
take as free parameters that we determine from our experimental binding energies. We perform a 

minimization with respect to the free parameters, and this allows us to determine the positions 
of both Feshbach resonances, and the direct mapping of the scattering length on the magnetic field 
a{B). In order to extract the resonance's parameters we fit a{B) with a factorized expression js?]: 
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where o^g is the background scattering length, Aj is the i's resonance width and i?o,j is the i's 
resonance position. 

We analyze the data in two different ways. First we perform an independent fit of the model to 
the Imp = 0) and the Imp = 1) binding energies, represented by solid blue and red lines in Fig.[71^a) 
and (b), respectively. Then, a combined fit for both mp states is performed (dashed black line in 
Fig. [71^a,b)), which allows us to check the consistency between the two different experiments, within 
the same interaction model. Deviations of the second case from the first one are a hint for possible 
different systematic shifts in magnetic field calibration for both states. This can be explained by the 
fact that a precise calibration of the magnetic field is performed locally by microwave transitions in 
the vicinity of each of the wide resonances. Note that each set of measurements is limited to a range 
of ±20 G around the resonances while the distance between them is ~ 150 G. We therefore use the 
results of the independent fits for the analysis which was discussed previously in this paper. 
The resonances' parameters as obtained from the independent fits are summarized in Table ITIl 



state 




type 


a-bg/ao 


A(G) 


Bo (G) 


Imp = 


0) 


narrow 


-18.24 


+4.518 


846.0 


Imp = 


0) 


wide 


-18.24 


-237.8 


893.7(3) 


\mp = 


1) 


wide 


-20.98 


-171.0 


738.2(2) 



Table II: Feshbach resonance parameters for both states as obtained from independent fits of the CC calcu- 
lation to the molecular binding energies (solid lines in Fig. [?])• 



Table IIIII represents the resonances' positions obtained from the combined fit and compares 
them with our most precise to-date experimental values, i.e. an atom loss measurement for the 
narrow resonance of \mp = 0) [6] and the two independent fits of the wide resonances (Table HI]). 
A field calibration uncertainty of 0.3 G is added to the fitting errors of the experimental values. 
We find the results in Table |TTI] mutually consistent. Moreover, the determined values for 4>s, 



A(j)T and Cg are all consistent with the bound of earlier performed analysis [18|, |29| . 



36|. Since 



our combined analysis only weakly depends on the value of the Cg coefficient, we present our results 
with a fixed value of Cq which was found in the ab initio calculations jsol . [sil . The determined 
parameters (l)s and Acpp correspond to the singlet and triplet scattering lengths 05=34.33(2) qq 
and aT=-26.87(8) qq. 
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state 




type 


Bo 

Combined fit 


(G) 

Experimental 


|tof = 


0) 


narrow 


845.54 


844.9(8) 


Imp = 


0) 


wide 


893.95(5) 


893.7(4) 


Imp = 


1) 


wide 


737.88(2) 


738.2(4) 



Table III: Feshbach resonances' positions (Bq) as obtained from a combined fitting of the molecular binding 
energy measurements in both states simultaneously to the CC calculation. The experimentally determined 
positions are presented in the last column where the narrow resonance was determined by atom loss mea- 
surement [sl and the two wide ones where determined by independent fits of the molecular binding energies 
(see Table HI]). The field calibration uncertainty is included in the experimental errors. 



IV. CONCLUSIONS 

In this work we study experimentally Efimov physics in two nuclear-spin sublevels of bosonic 
lithium and show that the positions and widths of recombination minima and Efimov resonances 
are identical for both states within the experimental errors. As the properties of Efimov features 
are governed by the three-body parameter, our study indicates that the short-range physics is 
nuclear-spin independent. We also find that the Efimov features are universally related across 
the Feshbach resonances. We note that slight deviations of our measurements from the universally 
predicted values can be explained by finite effective range corrections which were recently evaluated 



by means of an effective field theory 3a]. Let us note also that the observed position of the Efimov 
resonance reveals the same numerical factor |a_[/ro ~ 8.5 as in the experiments on ^'^'^Cs [2] which 
may or may not be an accidental coincidence. 

The reported results crucially depend on a careful mapping between the scattering length and 
the magnetic field. We characterize two wide Feshbach resonances in different states by fitting the 
binding energies of weakly bound molecules, created by rf-association, with a coupled channels 
analysis. This gives rise to a very precise determination of the absolute positions of the Feshbach 
resonances and the values of the singlet and triplet scattering length that characterize the molecular 
potentials of lithium. 
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